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Abstract 

We consider the problem of describing the asymptotic behaviour of FRW uni- 
KA , verses near their spacetime singularities in general relativity. We find that the 

j^ ■ Bel-Robinson energy of these universes in conjunction with the Hubble expansion 

rate and the scale factor proves to be an appropriate measure leading to a com- 
plete classification of the possible singularities. We show how our scheme covers 
all known cases of cosmological asymptotics possible in these universes and also 
predicts new and distinct types of singularities. We further prove that various 
asymptotic forms met in flat cosmologies continue to hold true in their curved 
counterparts. These include phantom universes with their recently discovered big 
rips, sudden singularities as well as others belonging to graduated inflationary 
models. 



1 Introduction 

The problem of cosmological singularities is perhaps the single most important unsolved 
issue in modern theoretical general relativity and cosmology. It is usually formulated 
as being a two-faceted problem, the first aspect being the issue of finding conditions 
under which singularities are expected to develop in a finite time (existence) during the 
evolution of a model spacetime, the second being that of describing their character or 
nature. In fact, the other, closely related, issue in general relativity, that of proving 
geodesic completeness, spacetimes existing for an infinite proper time interval, is usually 
thought of as the negation of the singularity problem, cf. [T. The geodesic completeness 
problem may also be formulated as a two-sided issue, one aspect being, like the singularity 
problem, that of finding conditions that guarantee the existence of complete spacetimes, 
the other being that of unraveling the nature of geodesic completeness, that is the nature 
of behaviour of the controlling functions, for instance, their decay to infinity. 

The issue of proving existence of singular and/or complete spacetimes is basically 
a geometric one. Traditionally one is interested in formulating criteria, of a geometric 
nature, to test under what circumstances singularities, in the form of geodesic incom- 
pleteness, will be formed during the evolution in general relativity and in other metric 
theories of gravity, cf. [2j. These, apart from the causality ones like global hyperbolicity, 
are usually translated into sufficient conditions to be satisfied by the matter fields present 
and are very plausible. Such conditions guarantee that in a finite time in the future or 
past of the spacetime, some quantity, usually describing the convergence of a bundle of 
geodesies, blows up thus designating the past or future breakdown of the bundle. On the 
other hand, one can formulate equally plausible and generic geometric criteria for the 
long time existence of geodesically complete, generic spacetimes in general relativity and 
other theories of gravity, cf. j^. Such sufficient conditions assume a globally hyperbolic 
spacetime in the so-called generic sliced form (cf. Q), require the space gradient of the 
lapse function as well as the extrinsic curvature not to grow without bound, and are also 
very plausible. They provide us with assumptions under which generic spacetimes will 



exist forever. All known global-in-time spacetimes satisfy these sufficient conditions and 
are therefore future geodesically complete. 

Of course as the singularity theorems do not detail the nature of the behaviour 
of the controlling functions on approach to the singularity, similarly the completeness 
theorems do not provide conditions for, say, the decay of such functions to infinity. A 
lesson to be drawn from this state of affairs is that in general relativity (and in other 
metric gravity theories), singular and complete spacetimes are equally generic in a sense. 
The real question is what do we mean when we say that a relativistic model develops 
a spacetime singularity during its evolution. In other words, what are the possible 
spacetime singularities which are allowed in gravity theories? These may be singularities 
in the form of geodesic incompleteness developing in the course of the evolution but 
can be also others which are more subtle and which are also dynamical and present 
themselves, for instance, in some higher derivatives of the metric functions spoiling the 
smoothness of global solutions to the field equations. 

Obviously a recognition and complete analysis of such a program cannot be accom- 
plished in the short run and requires complete examination of a number of different 
factors controlling the resulting behaviour. For instance, one needs to have control of 
each possible behaviour of the different families of relativistic geometry coupled to matter 
fields in general relativity and other theories of gravity to charter the possible singular- 
ity formation. In this sense, unraveling the nature and kinds of possible singularities in 
the simplest kinds of geometry becomes equally important as examining this problem 
in the most general solutions to the Einstein equations. In fact, an examination of the 
literature reveals that the more general the spacetime geometry considered (and thus 
the more complex the system of equations to be examined) the simplest is the type of 
singularities allowed to be examined. By starting with a simple cosmological spacetime 
we allow all possible types of singularities to come to surface and be analyzed. 

In |S] we derived necessary conditions for the existence of finite time singularities in 
globally and regularly hyperbolic isotropic universes, and provided first evidence for their 



nature based entirely on the behaviour of the Hubble parameter (extrinsic curvature). 
This analysis exploited directly the evidence provided by the completeness theorems 
proved in j3]. The main result proved in |3] may be summarized as follows: 

Theorem 1.1 Necessary conditions for the existence of finite time singularities in glob- 
ally hyperbolic, regularly hyperbolic FRW universes are: 

Si For each finite t, H is non-integrable on [ti,t], or 

S2 H blows up in a finite future time, or 

5*3 H is defined and integrable for only a finite proper time interval. 

Condition Si may describe different types of singularities. For instance, it describes 
a big bang type of singularity, when H blows up at ti since then it is not integrable on 
any interval of the form [ti,t], t > ti (regular hyperbolicity is violated in this case, but 
the scale factor is bounded from above). However, under Si we can have other types of 
singularities: Since H{t) is integrable on an interval [ti,t], if H{t) is defined on [ti,t], 
continuous on (ti, t) and the limits lim^^j+ H{t) and lim^-^i- H{t) exist, the violation of 
any of these conditions leads to a singularity that is not of the big bang type discussed 
previously. 

Condition 5*2 describes a future blow up singularity and condition 5*3 may lead to 
a sudden singularity (where H remains finite), but for this to be a genuine type of 
singularity, in the sense of geodesic incompleteness, one needs to demonstrate that the 
metric is non-extendible to a larger interval. 

Note that these three conditions are not overlapping, for example 5*1 is not implied 
by ^2 for if H blows up at some finite time tg after ti, then it may still be integrable on 
[ti,t], ti < t < tg. As discussed in [^, Theorem fll.l|) describes possible time singularities 
that are met in FRW universes having a Hubble parameter that behaves like Si or 5*2 
or 5*3. Although such a classification is a first step to clearly distinguish between the 
various types of singularities that can occur in such universes, it does not bring out some 



of the essential features of the dynamics that differ from singularity to singularity. For 
instance, condition S2 includes both a collapse singularity, where the scale factor a -^ 
as t — *> ts, and a blow up singularity where a ^ 00 as t ^ ts- Such a degeneration is 
unwanted in any classification of the singularities that can occur in the model universes 
in question. 

It is therefore necessary to extend and refine this classification by considering also 
the behavior of the scale factor. That it is only necessary to include this behaviour, is 
also seen most clearly by noticing that only in this way one may consistently distinguish 
between initial and final singularities, or between past and future ones, as it has been 
repeatedly emphasized by Penrose over the years in his Weyl curvature hypothesis, cf. [Hj. 
Another aspect of the problem that must be taken into account is the relative behaviour 
of the various matter components as we approach the time singularities. We believe that 
a clear picture of the behaviour of the matter fields on approach to a singularity will 
contribute to an understanding of its nature. 

In this paper we present a new scheme that describes in a consistent way the possible 
spacetime singularities which can form during the evolution of the FRW cosmologi- 
cal models in general relativity. Section 2 introduces the classification and shows that 
by extending our previous results and taking into account the effects induced by the 
Bel-Robinson energy, we can arrive at the possible forms of singular behaviour that 
the universes in question can exhibit. Section 3 proves that various singularity types 
predicted by our scheme have a general significance in the isotropic category and also 
describes how some singularity types known to occur only in flat isotropic universes do 
in fact appear in many curved spacetimes. In the last Section we discuss various aspects 
of this work including the issue of how to describe known types of singularities with our 
technique. 



2 Classification 

In this Section we present a classification of possible singularities which occur in isotropic 
model universes. Our classification is based on the introduction of an invariant geometric 
quantity, the Bel- Robinson energy (cf. and Refs. therein), which takes into account 
precisely those features of the problem, related to the matter contribution, in which 
models still differ near the time singularity while having similar behaviours of a and H. 
In this way, we arrive at a complete classification of the possible cosmological singularities 
in the isotropic case. 

The Bel-Robinson energy is a kind of energy of the gravitational field projected in a 
sense to a slice in spacetime. It is used in |H], [0] to prove global existence results in the 
case of an asymptotically flat and cosmological spacetimes respectively, and is defined 
as follows. Consider a sliced spacetime with metric 

{n+i)g = _^2^Q0Y^g,,QiQj^ ^o = t/t, e' = dx' + f3'dt, (2.1) 

where N = A^(t, x*) is the lapse function and (3'^{t,x^) is the shift function, and the 
2-covariant spatial electric and magnetic tensors 
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where rjij^ is the volume element of the space metric g. These four time-dependent space 
tensors comprise what is called a Bianchi field, {E,H,D,B), a very important frame 
field used to prove global in time results, cf. [7j. 
The Bel-Robinson energy at time t is given by 



B{t) 



I- [ {\E\' + \D\' + \B\' + \H\')di,-,„ (2.2) 

-^ JMt 



where by |Xp = g^^g^^X-ikXji we denote the spatial norm of the 2-covariant tensor X. 
In the following, we exclusively use an FRW universe filled with various forms of matter 
with metric given by 

ds^ = -dt^ + a^{t)da'^, (2.3) 

where da'^ denotes the usual time-independent metric on the 3-slices of constant curvature 
k. For this spacetime, we find that the norms of the magnetic parts, \H\^ \B\, are 
identically zero while \E\ and \D\, the norms of the electric parts, reduce to the forms 

|E|2 = 3(a/a)2 and \D\^ = 2> {{d/af + k/a^f . (2.4) 

Therefore the Bel-Robinson energy becomes 

B{t) = ^{\E\' + \D\'), (2.5) 

where C is the constant volume of (or in in the case of a non-compact space) the 3- 
dimensional slice at time t. 

It is not difficult to show that a closed, RW universe such that |-D| is bounded above, 
H must be bounded above and the scale factor bounded below. Therefore H must be 
integrable and the spatial metric bounded below, that is such a universe is regularly 
hyperbolic. This in turn means that all hypotheses of the completeness theorem proved 
in jH] are satisfied and therefore such a universe is f?— complete. It is also straightforward 
to see that the null energy condition is equivalent to the inequality \E\ < \D\, hence 
completeness is then accompanied with \E\ being bounded above. (For a fiat spatial 
metric we would need to impose the regular hyperbolicity hypothesis in order to conclude 
completeness since the latter is independent from the boundedness of \D\.) 

We can now proceed to list the possible types of singularities that are formed in an 
FRW geometry during its cosmic evolution and enumerate the possible types that result 
from the different combinations of the three main functions in the problem, namely, the 
scale factor a, the Hubble expansion rate H and the Bel Robinson energy B. These types 
wiU by necessity entail a possible blow up in the functions \E\, \D\. If we suppose that 



the model has a finite time singularity ai t = tg, then the possible behaviours of the 
functions in the triplet {H, a,{\E\,\D\)) in accordance with Theorem ()1.1|) are as follows: 

Si H non-integrable on [ti,t] for every t > ti 

5*2 H ^ oo at ts > h 

S3 H otherwise pathological 

Ni a^O 

N2 a —>■ as y^ 

iVs a — * 00 

Bi \E\ -^ 00, \D\ -^ cxo 
B2 \E\ < cx), \D\ -^ 00 
i?3 l^l -^ CX), \D\ < 00 
B4 \E\ < CX), \D\ < CX). 

The nature of a prescribed singularity is thus described completely by specifying the 
components in a triplet of the form 

with the indices i,j, I taking their respective values as above. 

Note that there are a few types that cannot occur. For instance, we cannot have an 
(5*2, N2, -B3) singularity because that would imply having a < cx (A''2) and H ^ 00 {S2), 
while 3 ((d/a)^ + k/a^) < cx (-B3), at tg which is impossible since |Z)p ^ cx at t^ [k 
arbitrary) . 

A complete list of impossible singularities is model-dependent and is generally given 
by triplets {Si, Nj, Bi), where the indices, in the case of a A; = 0, +1 universe, take the 



values i = 1,2, j = 1, 2, 3, / = 3, 4, whereas in the case of a A; = —1 universe the indices 
take the values i = 1,2, j = 2,3, I = 3,4 (here by 5*1 we denote for simplicity only 
the big bang case in the 5*1 category). We thus see that some singularities which are 
impossible for a flat or a closed universe become possible for an open universe. Consider 
for example the triplet (5*2, A^^i, B^) which means having H -^ oo, a ^ and 



\D\^ = 3{id/af + k/a^Y< 



CXD 



at ts- This behaviour is valid only for some cases of an open universe. 

All other types of flnite time singularities can in principle be formed during the 
evolution of FRW, matter-flUed models, in general relativity or other metric theories of 
gravity. 

It is interesting to note that all the standard dust or radiation-fllled big bang sin- 
gularities fall under the strongest singularity type, namely, the type (5*1, A^^i, -Bi). For 
example, in a flat universe fllled with dust, at t = we have 

a{t) oc t2/3 -> 0, (iVi), (2.6) 

H oc t-^ ^ oo, (Si), (2.7) 

|E|2 = 3/4H*^oo, IDI"^ = 3H^ ^ oo, (Bi). (2.8) 

Note that our scheme is organized in such a way that the character of the singularities 
(i.e., the behaviour of the deflning functions) becomes milder as the indices of S, N and 
B increase. Milder singularities in isotropic universes are thus expected to occur as one 
proceeds down the singularity list. 

It is the purpose of this classification to apply both to vacuum as well as matter 
dominated models. In fact the Bel-Robinson energy takes care in a very neat way the 
matter case. For instance, in fluid-fllled models, the various behaviours of the Bel- 
Robinson energy density can be related to four conditions imposed on the density and 
pressure of the cosmological fluid: 

Bi <^ /i — !► oo and |/i + 3p| — *> oo 



B2 -^ fi -^ oc and |/i + 3^1 < oo 

-B3 <^ /i < 00 and l/i + 3p| -^ 00 <^ /i < CX3 and \p\ -^ 00 

i?4 -vv- /x < 00 and l/x + 3p| < 00 -vv- /I < 00 and |p| < 00. 

Of course we can translate these conditions to asymptotic behaviours in terms of a, H, 
depending on the value of k, for example, 

1. If fc = 0, yU < 00 ^ H'^ < 00, a arbitrary 

2. If /c = 1, yU < 00 ^ i/^ < 00 and a 7^ 

3. If A; = -1, n < 00 ^ H'^ - l/a^ < cx). 

As an example, we consider the sudden pressure singularity introduced by Barrow in 
[TU] . This has a finite a (condition N2), finite H (condition 5*3), finite /z but a divergent 
p (condition B3) at tg. 

As another example, consider the flat FRW model containing dust and a scalar field 
studied in jTT]. The scale factor collapses at both an initial (big bang) and a final time 
(big crunch). The Hubble parameter and a/ a both blow up at the times of the big bang 
and big crunch (cf. f5|) leading to an {Si, Ni, Bi) big bang singularity and an (S'2, Ni, Bi) 
big crunch singularity, respectively. 

3 Generic results and examples 

It was recently shown by Ellis in ^21 that a RW space with scale factor a{t) admits a 
past closed trapped surface if the following condition is satisfied: 



a{t) > 



(3.1) 



fir) 

with /(r) = sinr, r, sinhr for k = 1,0,-1 respectively. Recall that a closed trapped 
surface is a 2-surface with spherical topology such that both families of incoming and 



outgoing null geodesies orthogonal to the surface converge. Since our function \D\ can 
be written in the form 



IDl ^ 



a2(t) 



m-m)(m'^'^'^ 



(3.2) 



fir) A fir) J f\r) 

we see that the condition for the existence of a closed trapped surface becomes equivalent 
to the following inequality: 

We thus conclude that collapse singularities (as predicted by the existence of a trapped 
surface) are characterized by a divergent Bel-Robinson energy. 

In this Section we provide necessary and sufficient conditions for the occurrence of 
some of the triplets detailing the nature of the singularities introduced above. These 
conditions are motivated from studies of cosmological models described by exact solu- 
tions in the recent literature. By exact solutions we mean those in which all arbitrary 
constants have been given fixed values. We expect the proofs of these results to be all 
quite straightforward, for we have now already identified the type of singularity that we 
are looking for in accordance with our classification. Proving such results without this 
knowledge would have been a problem of quite a different order. 

The usefulness of the results proved below lies in that they answer the question of 
whether the behaviours met in known models described by exact solutions (which as a 
rule have a flat spatial metric {k = 0)) continue to be valid in universes having nonzero 
values of k as well as described by solutions which are more general than exact in the 
sense that some or all of the arbitrary constants present remain arbitrary. The reason 
behind this behaviour whenever it is met, lies in the fact that the curvature term in the 
Friedman equation turns out to be usually subdominant compared to the density term 
or in any case can not alter the H behaviour. For the purpose of organization we present 
separately results about future and past singularities. 
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3.1 Future singularities 

The first result of tliis subsection cliaracterizes tlie future singularity in phantom cos- 
mologies irrespective of the value of the curvature k and says that the singularities in 
such models can be milder than the standard big crunches and have necessarily diverging 
pressure and the characteristic "phantom" equation of state. 

Theorem 3.1 Necessary and sufficient conditions for an {S2,N-^,Bi) singularity occur- 
ring at the finite future time tg in an isotropic universe filled with a fluid with equation 
of state p = wfi, are that w < —1 and \p\ ^ oo at tg. 

Proof. Substituting the equation of state p = w^ in the continuity equation /i + 
3H{fi + p) = 0, we have 

/i oc a-^^'^+^l (3.4) 

and so if w < — 1 and p blows up at t^, a also blows up at tg. Since 

H' = ^-^,^ \D\' = Y^ |Er = lAl + 3u;)^ (3.5) 

we conclude that at tg, H, a, \D\ and \E\ are divergent. 

Conversely, assuming an (5*2, A^s, -Bi) singularity at tg in an FRW universe with the 
equation of state p = wfi, we have from the [Bi) hypothesis that fi -^ oo a.t tg and so p 
also blows up at tg. Since a is divergent as well, we see from ()3.4|) that w < —1. m 

As an example, consider an exact solution which describes a flat, isotropic phantom 
dark energy filled universe, studied in |12|| given by 



a = [ao'^'-^-y^ + ^ii±|M(t - to)]''""' , (3.6) 

where A is a constant. From ()3.4j) we see that the scale factor, and consequently H, 
blows up at the finite time 

2 

^s — ^0 + 



3VA{\w\ - l)ao 
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3(|«;|-l)/2- 



Then it follows that |Ep = |i7^(l + Sw)^ and \D\'^ = 3H* also blow up at t,. Therefore 
in this model the finite time singularity is of type (5*2, A^s, -Bi). 

Next we focus on a generalization of a model, called graduated inflation and originally 
proposed by Barrow - see next subsection, given in f^^. Consider a fiat FRW filled with 
a fluid with equation of state p + /i = —B^^, [3 > 1. As /x ^ oo, 

t -^ to + 



V3k5 1-2/3 

where k^ = 8ttG and to is an integration constant. It follows then that \p\ -^ oo at 
to- The scale factor is described by ()3.7p below and it is therefore finite. However, 
H, \D\ and \E\ all diverge at to meaning that this is an {S2, N2, Bi) singularity in our 
scheme. As we shall show in the next subsection, the introduction of a more general 
equation of state has resulted in "taming" the singularity of the graduated inflationary 
type. This behavior is a special case of the following general result which holds also in 
curved models. 

Theorem 3.2 A necessary and sufficient condition for an (5*2, A^2, -Bi) singularity atts 
in an isotropic universe filled with a fluid with equation of state p + ft = —B^i^, [3 > 1, 
is that /i ^ 00 attg. 



Proof. From the continuity equation we have 

.3^(1-/3)7' 
and so a ^ Oq as t — > tg. Since 



a = ao exp ( — — — ) , (3.7) 



^^ = f-^, |i?P = l(2/. + 3i?/)^ \D? = ^, (3.8) 

we see that as t ^ t^, if, \E\ and \D\ diverge, leading precisely to an (5*2, A^2, -Bi) 
singularity. The converse is obvious. ■ 
Similarly one can prove: 
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Theorem 3.3 A necessary and sufficient condition for an (S'3, A^2, -S3) singularity atts 
in an isotropic universe filled with a fluid with equation of state p + fJ, = —C{fj,Q — fj.)~'^ , 
J > 0, is that fi ^ fiQ at ts- 



aocexp<'-^^^^— ^|>, (3.9) 



Proof. Again using the continuity equation we find 

(/JQ - /i)^+ 

3C(7 + 1) 
which is finite as t ^ t^. Also since 

^^ = f-^, |i?P = l(2/. + 3C(/.o-/i)--)^ \D\' = ^, (3.10) 

we see that as t —>■ tg, H and \D\ remain finite whereas \E\ diverges, leading to an 
{S3, N2, B^i) singularity. The converse is immediate. ■ 

The above equation of state is studied in jHj for the case of a fiat universe. It follows 



then that as ;U — > /io, t = to — ^^cva" ( +1) ~^ ^^ ^^^ integration constant) and \p\ -^ 00. 
The scale factor is described by (|3.9|) and it is therefore finite at to whereas H < 00, 
\D\ < 00 and \E\ — > 00 leading to an (6*3, A^2! -B3) singularity. We therefore see that 
the nature of the singularity depends very sensitively on even very mild changes in the 
equation of state. 

A way to become more intimately acquainted with the nature of the various singulari- 
ties and to clearly distinguish their various differences is to study the relative asymptotic 
behaviours of the three functions that define the type of singularity. Using a standard 
notation that expresses the behaviour of two functions around the singularity at t*, we 
can introduce a kind of relative "strength" in the singularity classification. Let /, g be 
two functions. We say that 

1. f{t) is much smaller than g{t), f(t) << g(t), if and only if 
limt^t, f{t)/g{t) = 

2. fit) is similar to git), fit) ~ git), if and only if < lim^^t^ fit) /git) < 00 

3. fit) is asymptotic to git), fit) <-> git), if and only if limi^j, fit) /git) = 1, 
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Using these notions we find that standard radiation filled isotropic universes (with 
k = 0, ±1) have the asymptotic behaviours described by a << H « {\E\ *-^ |-D|), 
whereas the rest of the standard big bang singularities have a « H « {\E\ ~ \D\). 

The phantom model of theorem (jH.lj) has three possible behaviours depending on the 
ranges of the w parameter, namely, if —5/3 < w < —4/3 we have H « a « {\E\ ~ 
\D\), if -4/3 < w < -I then H « {\E\ ~ \D\) « a whereas if u; < -5/3 we have 
a«H « {\E\ ~ \D\). 

The fiuid case of theorem ()3.2|1 has the behaviour is a << H << \D\ << \E\. 

The sudden singularity of theorem ()3.3|) is characterized by the behaviour {H ~ |Z}| ~ 
a) << \E\. If we identify a sudden singularity as the kind of singularity occurring with 
a divergent p while a and H remain finite, then this asymptotic behaviour is the only 
possibility. 

In the majority of the types of singularities that we have met in our study the two 
quantities \E\, \D\ are the most wildly diverging functions and therefore become the 
most dominant ones asymptotically. 

3.2 Past singularities 

We now present two results on the nature of past time singularities. These include the 
standard big bang ones but the latter are obviously not the only possibility. The first 
result given below predicts that the time singularity met in an exact solution in the 
family of the so-called graduated infiationary models first constructed in ^5], extends 
also to open universes. The model consists of a fiat isotropic spacetime with a fiuid with 
equation of state p + /i = 7/^^^"^, 7 > 0, and admits the exact solution 

/^ 16 \ ., 256 ,^^^, 

« = ^^P (^-^375:^J ' ^ith /^=^:^- (3.11) 

At t = 0, a -^ (Ni), H ^ oo (Si) and jEp = ^(-2/i + 3-ffj,^/y -^ oo, jDp -^ oo 
(Bi). This behaviour can be made to occur more generally due to the following result. 
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Theorem 3.4 A necessary and sufficient condition for an {Si,Ni,Bi) singularity at ti 
in an open or fiat isotropic universe filled with a fluid with equation of state p + fj, = ■jfj,^, 
7 > and X < 1, is that fj, ^ oo at ti. 



Proof. The continuity equation gives 



(j^'^' 



» = "°"^P(^ 37(A-1) J' ^^-^^^ 



so that a ^ as t ^ ti. Since 






|i^r = ^(-2^ + 37/i')', 



II 3 , 



we see that as t — > ti H, \E\ and \D\ diverge provided that fc = or A; = —1. The 
converse is straightforward. ■ 

The asymptotic strength of this singularity is a << {\E\ ^^ \D\) << H . 

The second resuh that we prove now says that the strongest big bang type singular- 
ities in universes with a massless scalar field are produced due to the special form of the 
kinetic term. An exact solution for the case of a flat spatial metric is derived in ^j and 
is given by 

3t' ^ V 3 c 

Since a oc t^/^ we have that at t = 0, a -^ (A^i), H = l/(3t) -^ 00 (Si), \E\^ = 073 -* 
00 and |i5p = 3H^ -^ 00 (-Bi). As it follows from fl^'l, this exact solution represents the 
asymptotic behaviour of a scalar field model if 

lim e-^l'^l\/(0) = 0. 

fji— >itoo 

Theorem 3.5 A necessary and sufficient condition for an {Si, Ni, Bi) singularity at ti 
in an isotropic universe with a massless scalar field is that — > cx) at ti. 
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Proof. From the continuity equation, (f) + 3H(j) = 0, we have oc a ^, and if -^ oo 
then a ^ 0. Since 

H' = ^--,^oo (3.13) 

H becomes unbounded at ti. In addition, since 

2^/?^0^^ (3.14) 

' ' 3 12 ' ^ ' 



and 



|Ep = ^(/i + 3p)2 = ^^oo, (3.15) 



at ti, both \D\ and \E\ diverge there. 

Conversely, assuming an (5*1, Ni,Bi) singularity at ti, we have (from Bi) that /x ^ oo 
and so 0^ ^ oo at ti. ■ 

The asymptotic strength of the singularity for the model described by both the exact 
solution and the above theorem is a << H << {\E\ ~ |-D|). 

A very mild type of singularity discovered recently in |T3]. We shall see how our 
classification is able to pick this up. The model is that of a flat FRW type filled with a 
fluid with equation of state 

and for /? = 1/5 it admits the exact solution a = a^e^ . Then H = {8/3)t^''^, H = 
(40/9)r^/^ and H = (80/27)r~^/^. Thus as r -^ 0, a, d, d, H, H all remain finite whereas 
H becomes divergent. 

We can easily see that in this universe the Bel Robinson energy at the initial time, 
B{0), is finite whereas its time derivative is 



i3(. 



2^(^ + 2M)+4(^ + i.^)(-^ + ™) 



(3.16) 



and thus B{t) ^ oo at r — i> 0. This is like a B4 singularity in our scheme. Since 
the derivative of the Bel-Robinson energy diverges, we may interpret this singularity 
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geometrically as one in the velocities of the Bianchi (frame) field. At tg the Bianclii field 
encounters a cusp and its velocity diverges there. We believe this to be the mildest type 
of singularity known to date. 

4 Discussion 

In this paper we have extended and refined the classification of the cosmological singu- 
larities possible in an FRW universe in general relativity. A classification based entirely 
on the asymptotic behaviour of the scale factor or the Hubble expansion rate cannot 
lead to complete results and we have found it necessary to extend this scheme to include 
the behaviour of the Bel-Robinson energy. In this case we have found that the result- 
ing behaviours of the three functions, H, a and B taken together exhaust the types of 
singularities that are possible to form during the evolution of an isotropic universe. 

The resulting classification is described by triples of the form {Si,Nj,Bi) where the 
S category monitors the asymptotic behaviour of the expansion rate, closely related to 
the extrinsic curvature of the spatial slices, A^ that of the scale factor, describing in a 
sense what the whole of space eventually does, while B describes how the matter fields 
contribute to the evolution of the geometry on approach to the singularity. We know (cf. 
[SlITj) that all these quantities need to be uniformly bounded to produce geodesically 
complete universes. Otherwise, the whole situation can be very complicated and we have 
exploited what can happen in such a case when we consider a relatively simple geometry 
as we have done in this paper. 

Our scheme not only does cover all the recently discovered types of singularities but it 
also predicts many possible new ones. For example, in the case of a fiat isotropic universe 
the classification of [14J provides us with four main types of singularity. These four types 
can be identified with four particular {Si,Nj,Bk) triplets of our scheme. The 'big rip' 
type characterized by a — > oo, /x — > oo and \p\ ^ oo at t^ is an (S'2, A3, Bi) singularity; 
the 'sudden' singularity described by a — »• a^ < 00, /x — >• /i^ < cxd and \p\ — > 00 at t^ is 
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an {S^, N2, B3); further, the so-called type III singularity, namely a ^ a.^ < 00, ;U — » 00 
and \p\ ^ 00 at ts is clearly an (5*2, N2, Bi) singularity, while, type IV singularities with 
a ^ a^ < 00, /i — * and \p\ ^ at t^ belong to an (5*3, N2, -B4). 

Having all possible singularity types in the form expounded in this paper has the 
added advantage that we can consistently compare between different types as we asymp- 
totically approach the time singularity. In this case the relative strength of the functions 
describing the singularity type becomes an important factor to distinguish between the 
possible behaviours. In this way we can have a clear picture of how "strong" a singu- 
larity can be when formed during the evolution and also why such a type eventually 
arises in any particular model. This charters the singularities in the isotropic category, 
a necessary first step in an attempt to consider the same classification problem in more 
complex situations. 

It is natural to consider the extension of the work done in this paper in the context 
of the more general anisotropic Bianchi models. We believe that an analysis of this 
more complicated case is still feasible using the techniques of the present paper. In such 
contexts the topology and spatial extent of the singular surfaces is expected to play a role 
even for the simplest non-trivial 'vacua' such as the Kasner or the Taub-NUT solutions 
and a complete elucidation of these cases is a prerequisite to situations with matter fields. 
We leave such matters to future work. 
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